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We present strange particle spectra and yields measured at mid-rapidity in ^Js = 200 GeV proton- 
proton (p + p) collisions at RHIC. We find that the previously observed universal transverse mass 
(mx = \J pt 2 + m 2 ) scaling of hadron production in p + p collisions seems to break down at higher 
hit and that there is a difference in the shape of the hit spectrum between baryons and mesons. 
We observe mid-rapidity anti-baryon to baryon ratios near unity for A and E baryons and no 
dependence of the ratio on transverse momentum, indicating that our data do not yet reach the 
quark-jet dominated region. We show the dependence of the mean transverse momentum ((pt)) on 
measured charged particle multiplicity and on particle mass and infer that these trends are consistent 
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with gluon-jet dominated particle production. The data are compared to previous measurements 
from CERN-SPS, ISR and FNAL experiments and to Leading Order (LO) and Next to Leading 
order (NLO) string fragmentation model predictions. We infer from these comparisons that the 
spectral shapes and particle yields from p + p collisions at RHIC energies have large contributions 
from gluon jets rather than quark jets. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.40.Ep 



I. INTRODUCTION 

The production of particles in elementary proton- 
proton (p+p) collisions is thought to be governed by two 
mechanisms. Namely, soft, thermal-like processes which 
populate the low momentum part of the particle spectra 
(the so-called underlying event) and the hard parton- 
parton interaction process. In this scenario, the low 
transverse momentum (pt) part of the spectrum is expo- 
nential in transverse mass (mx= \J m 2 + pt 2 ) while frag- 
mentation, in leading order models, introduces a power 
law tail at high px. We investigate the validity of these 
assumptions at RHIC energies by studying the spectral 
shapes and the yields of identified strange hadron spec- 
tra from the lightest strange mesons (K^) to the heavy, 
triply-strange baryon. 

In this paper we report the results for transverse mo- 
mentum spectra and mid-rapidity yields (dN/dy) of , 
Kg, A, A, E-, S + , and fT +H + measured by the STAR 
experiment during the 2001-2002 -^8=200 GeV p+p run- 
ning at RHIC. After a brief introduction in Section UTI of 
the experimental setup and the conditions for this run, 
we provide a description of the event selection criteria 
and the efficiency of reconstructing the primary interac- 
tion vertex in Section IIII Al Specific attention will be 
given to the complications introduced by more than one 
event occurring in the detector during readout, a con- 
dition referred to as "pile-up". The details of strange 
particle reconstruction and the efficiency thereof will be 
discussed in Sections IIII Bl IIII CI and IIII Dl In Section 
II V Al we describe the final measured px spectra and mid- 
rapidity yields. We also describe the functions that were 
used to parameterize the px spectra in order to extrap- 
olate the measurement to zero px- We will show that 
the previously widely used power-law extrapolation for 
p + p and p + p collisions |l| does not yield the best 
X 2 results for the strange baryons and we will consider 
alternatives. Section FlV Bl introduces the idea of trans- 
verse mass scaling (mx-scaling) and its applicability to 
our data. The measured anti-particle to particle ratios 
are presented in Section TlV CI Interesting trends of in- 
creasing mean transverse momentum, (px), with particle 
mass have been previously observed in p + p collisions 
at ISR energies (20 < y/s < 63 GeV) Q Mean trans- 
verse momentum has also been found to increase with 
event multiplicity in p + p collisions at SppS (y/s = 630 
GeV) and FNAL energies (300 GeV < Vs < 1.8 TeV) 
0, 0. We will show the dependence of our (px) mea- 
surements on both particle mass and event multiplicity 



in Section IIVDI We discuss the details of the experi- 
mental errors and then compare our results in Section 
[3 with several models that attempt to describe particle 
production in p + p collisions via pQCD, string fragmen- 
tation, and mini-jets 0- We conclude in Section IVll with 
a discussion of the major results and some remarks about 
future directions for the ongoing analyses. 



II. EXPERIMENTAL SETUP 

The data presented in this paper were collected with 
the STAR detector [|J . The primary detector sub-system 
used for these analyses is the large cylindrical Time Pro- 
jection Chamber (TPC), which is able to track charged 
particles in the pseudo-rapidity range \rj\ < 1.8 with full 
azimuthal coverage 0. The TPC has 45 pad rows m 
the radial direction allowing a maximum of 45 hits to 
be located on a given charged particle track. A uniform 
magnetic field of 0.5 T is applied along the beam line 
by the surrounding solenoidal coils allowing the momen- 
tum of charged particles to be determined to within 2-7% 
depending on the transverse momentum of the particle. 
The field polarity was reversed once during the 2001- 
2002 run to allow for studies of systematic errors. The 
TPC tracking efficiency in p + p collisions is greater than 
90% for charged particles with Px> 300 MeV/c in the 
pseudorapidity region \rj\ < 0.7 0|. Particle identifica- 
tion may be achieved via measurements of energy loss 
due to specific ionization from charged particles passing 
through the TPC gas (dE/dx). The dE/dx, when plot- 
ted vs. rigidity separates the tracks into several bands 
which depend on the particle mass. A semi-empirical for- 
mula describing the variation of dE/dx with rigidity is 
provided by the Bethe-Bloch equation 0. An updated 
form, which accounts for the path length of a given par- 
ticle through matter, has been given by Bichsel and pro- 
vides a reasonable description of the dE/dx band centers 
for the particles presented in this paper 0. The Bichsel 
curves are shown in Figure ^ 

The dataset analyzed in this paper consisted of 1.4 x 
10 7 minimally-biased events before cuts. After applying 
a cut requiring the location of the primary vertex to be 
within 50 cm of the center of the TPC along the beam 
axis, to limit acceptance variations, 6 x 10 6 events re- 
mained. In all events, the detectors were triggered by re- 
quiring the simultaneous detection of at least one charged 
particle at forward rapidities (3.5 < \rj\ < 5.0) in Beam- 
Beam scintillating counters (BBCs) located at both ends 
of the TPC. This is referred to as a minimally biased 
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FIG. 1: dE/dx vs. Momentum for STAR p + p collisions at 
-^13=200 GeV. The curves are Bichsel parameterizations ||- 



trigger. The BBCs are sensitive only to the non-singly 
diffractive (NSD) part (30 mb) of the p+p total inelastic 
cross-section (42 mb) 0, QI|- A more detailed descrip- 
tion of STAR in general [6| and the complete details of 
the TPC in particular Jjj can be found elsewhere. 



III. ANALYSIS 
A. Primary Vertex Finding and Event Selection 

The position of the interaction vertex is calculated by 
considering only those tracks which can be matched to 
struck slats of the STAR central trigger barrel (CTB) 
[TT| . The CTB is a scintillating detector coarsely seg- 
mented into 240 slats placed azimuthally around the out- 
side of the STAR TPC at a radius of 2 m. It has a total 
pseudorapidity coverage of— 1.0 < r/ < 1.0 and has a fast 
response time of 10-60 ns, which is roughly one quar- 
ter of the time between beam bunch crossings (218 ns in 
the 2001-2002 run). Therefore, in approximately 95.6% 
of our p + p collisions, only charged particles from the 
triggered event will produce signals in the CTB which en- 
sures that the primary vertex is initiated with tracks from 
the triggered event only (note that, unlike the BBCs, the 
CTB itself is not used as a trigger detector for the event 
sample presented here). Furthermore, the primary ver- 
tex is assumed to be located somewhere along the known 
beam line. The z co-ordinate (along the beam) of the 
primary vertex is then determined by minimizing the \ 2 
of the distance of closest approach of the tracks to the 
primary vertex. 

The RHIC beams were tuned so as to maximize the lu- 
minosity and, consequently, the number of collisions that 
can be recorded. The average RHIC luminosities, which 
varied from 5 x 10 28 cm _2 s _1 to 5 x 10 30 cm _2 s~ 1 , pro- 



duce collisions more frequently (on the order of 2-200 
kHz) than the TPC can be read out (100 Hz). During 
p + p running, as many as five pile-up events can overlap 
(coming in the ^39 /is before or after an event trigger) 
in the volume of the TPC. Pile-up events come earlier 
or later than the event trigger and tracks from pile-up 
events may therefore be only partially reconstructed as 
track fragments. These track fragments from a pile-up 
event can distort the determination of the location of the 
primary interaction vertex as they do not point back to 
the vertex of the triggered event. To solve this, tracks 
that do not match to a struck CTB slat are not used in 
the determination of the primary vertex position. The 
remaining pile-up tracks, which match by chance to fired 
CTB slats, can then be removed with a reasonably re- 
strictive (2-3 cm) analysis cut on a track's distance of 
closest approach to the determined primary vertex. 

Another problem faced in the event reconstruction is 
the observation that for many minimally-biased triggers 
no primary vertex is reconstructed. The problem is sys- 
tematically worse for the low multiplicity events. There- 
fore, a correction must be applied to account for the 
events that are triggered on yet lost in the analyses due 
to an unreconstructed primary vertex. 

The efficiency of the primary vertex finding software 
was investigated by generating Monte Carlo (MC) p + p 
events, propagating the Monte Carlo produced particles 
through the STAR detector simulation (GEANT), then 
adding the resulting simulated signals into the abort-gap 
events. In an abort-gap event, the detectors are inten- 
tionally triggered when there are no protons in one or 
both of the beam bunches passing through the detec- 
tor. Abort-gap events therefore contain background due 
to the interaction of beam particles with remnant gas 
in the beam pipe and may also contain background re- 
maining in the TPC from collisions in the crossings of 
previous or subsequent beam bunches. Abort-gap events 
provide a realistic background environment in which to 
simulate the vertex finding process. The embedded simu- 
lated event is then passed through the full software chain 
and tracks are reconstructed. These events are then com- 
pared to the input from the MC events. A quantity A(z), 
representing the difference along the z (beam) axis be- 
tween the actual embedded MC primary vertex and the 
reconstructed primary vertex is defined as follows: 

A / \ LMC ^reconstructed I /i\ 

1\{Z) — Zpy — Zpy . (L) 



The probability distribution of A(z) is shown in Figure 
121 for approximately 87,000 simulated events. We sepa- 
rate events where the software finds a vertex into two 
classes. An event with a good primary vertex is defined 
as having A(z) < 2 cm, whereas a fake vertex event is 
one in which A(z) > 2 cm. Whilst this limiting value is 
somewhat arbitrary, it does relate to offline cuts in our 
particle reconstruction that are sensitive to the accuracy 
of the found vertex. 
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simulated events: 86733 
PV not reconstructed: 12160 
Good Vertex: A(z) < 2cm: 67014 
Fake Vertex: A(z) > 2cm: 7559 
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FIG. 2: Distribution of A(z). The unshaded region is the 
accepted range of good reconstructed event vetices. 



FIG. 3: Primary vertex finding efficiency vs measured primary 
track multiplicity. The horizontal line at one is only a guide 
for the eye. 



It was found that the probability of finding the pri- 
mary vertex was strongly dependent on multiplicity. For 
the purposes of this study, "charged track multiplicity" 
is defined as being a count of tracks in the TPC that 
have at least 15 hits, at least 10 of which must be used 
in the track fit. After separating the raw charged track 
multiplicity distributions for each event class, i.e. lost 
vertex, fake vertex and good vertex, these distributions 
can be divided by the charged track multiplicity distri- 
bution of all events. This ratio then represents the prob- 
ability for a certain event class to occur as a function of 
the measured charged track event multiplicity. Finally, 
the probabilities for each charged track multiplicity are 
mapped back to the corresponding primary track multi- 
plicity, where "primary tracks" are those which satisfy 
the above requirements and additionally point back to 
within 3 cm of the primary vertex. The probabilities for 
each event class as a function of primary track multiplic- 
ity are shown in Figure |3| Whereas lost vertex events 
are monotonically decreasing with increasing multiplic- 
ity, fake vertex events are most probable when the event 
has 2 primary tracks. The open symbols in Figure|3show 
the corresponding "found" and "good" probabilities for 
events that contained at least one strange particle decay 
candidate. Note that primary vertex finding is initiated 
with tracks pointing at fired slats of the CTB, as men- 
tioned above. But all found tracks are allowed to con- 
tribute to the final vertex position. Therefore, on rare 
occasions and in low multiplicity events, a vertex may be 
found with no single track pointing back within 3 cm. 
These events will appear in Figure [3] as having a found 
(or fake) vertex but zero primary track multiplicity. The 



use of these probabilities to correct the strange particle 
yields and event counts as a function of multiplicity is 
described later. 



B. Particle Identification 

All the strange particles presented here, with the ex- 
ception of the charged kaons, were identified from the 
topology of their weak decay products in the dominant 
channel: 

K° s 7T+ + 7T- (68.6%) (2) 



A -> p + n~ (63.9%) (3) 



3- -> A + 7T- (99.9%) (4) 



fi~ -> A + K~ (67.8%) (5) 

The charged tracks of the daughters of neutral strange 
particle decays form a characteristic "V" -shaped topolog- 
ical pattern known as a "VO" . The VO finding software 
pairs oppositely charged particle tracks to form VO candi- 
dates. These candidates can then be further paired with 
a single charged track, referred to as the "bachelor" to 
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form candidates for S~ and fl decays. During the ini- 
tial finding process, loose cuts are applied to partially 
reduce the background while maximizing the candidate 
pool. Once the candidate pool is assembled, a more strin- 
gent set of cuts is applied to maximize the signal-to-noise 
ratio and ensure the quality of the sample. The cuts are 
analysis dependent and are summarized in Table [I] for 
the Kg and A analyses and in Table ITT| for the S and f2 
analyses. 



Cut 


K° s 


A and A 


DC A of VO to primary vertex 


< 2.0 cm 


< 2.0 cm 


DCA of VO-daughters 


< 0.9 cm 


< 0.9 cm 


N(hits) daughters 


> 14 


> 14 


N(o-) dE/dx 


< 3 


< 5 


Radial Decay Length 


> 2.0 cm 


> 2.0 cm 


Parent Rapidity (y) 


± 0.5 


± 0.5 



TABLE I: Summary of Kg and A cuts. See text for further 
explanation. 



Cut 


H and H + 


tt and Cl + 


Hyperon Inv. Mass 


1321 ± 5 MeV 


1672 ± 5 MeV 


Daughter A Inv. Mass 


1115 ± 5 MeV 


1115 ± 5 MeV 


N(er) dE/dx bachelor 


< 5 


< 3 


N(<t) dE/dx pos. daugh. 


< 5 


< 3.5 


N(cr) dE/dx neg. daugh. 


< 5 


< 3.5 


N(hits) bachelor 


> 14 


> 14 


N(hits) pos. daugh. 


> 14 


> 14 


N(hits) neg. daugh. 


> 14 


> 14 


Parent Decay Length (lower) 


> 2.0 cm 


> 1.25 cm 


Parent Decay Length (upper) 


< 20 cm 


< 30 cm 


Daugh. V0 Decay Length (lower) 


N/A 


> 0.5 cm 


Daugh. V0 Decay Length (upper) 


N/A 


< 30 cm 


DCA of Parent to PVtx 


N/A 


< 1.2 cm 


DCA of Daughters 


N/A 


< 0.8 cm 


DCA of V0 Daughters 


N/A 


< 0.8 cm 


DCA of Bachelor to PVtx (lower) 


N/A 


> 0.5 cm 


DCA of Bachelor to PVtx (upper) 


N/A 


< 30 cm 


Parent Rapidity 


± 0.5 


± 0.5 



TABLE II: Summary of H and Q cuts. See text for further 
explanation. 



Several of these cuts require some further explanation. 
A correlation has been observed between the luminos- 
ity and the raw V0 multiplicity. This correlation is sug- 
gestive of pile-up events producing secondary VOs. The 
apparent path of the V0 parent particle (the Kg or A) 
is extrapolated back towards the primary vertex. The 
distance of closest approach (DCA) of the V0 parent to 
the primary vertex is then determined. Secondary VOs 



from pile-up events do not point back well to the pri- 
mary vertex of the triggered event and may therefore be 
removed via a cut on the DCA of the V0 parent to the 
primary vertex. Parent particles for secondary VOs may 
be charged and curve away from the primary vertex be- 
fore decaying, causing the secondary V0 to also point 
back poorly. Therefore, this cut also removes some true 
secondary VOs. 

Tracks in the TPC are occasionally broken into two or 
more segments that appear to be independent tracks to 
the V0 and S finding software. In the majority of cases, 
this is due to tracks crossing the boundaries between sec- 
tors of the TPC pad plane. A cut requiring a minimum 
number of hits is applied to each of the decay daughter 
tracks to minimize the contamination from these track 
fragments. 

Lastly, we define a variable N(cr) to quantitatively mea- 
sure the distance of a particular track to a certain particle 
band in dE/dx vs. rigidity space |12U13|: 



N(a) 



dE/dx 



measured 



dE/ (iXBichsel 



samples 



) • dE/dx 



(6) 



measured 



where R is the dE/dx resolution (width in dE/dx of the 
distribution of a given particle band, see Fig. ^) at the 
track's momentum and N samp i cs is the number of hits 
used in the determination of the dE/dx. Cutting on the 
N(cr) of a given track helps to decrease the background 
even further by decreasing the contamination of the can- 
didate pool due to misidentified tracks. This is particu- 
larly important for the Q analysis. The A and 5 analyses 
can tolerate more open cuts in favor of increased statis- 
tics. The invariant mass distributions for Kg, A, E — , il 
and their corresponding antiparticles are shown in Figure 

m 

The charged kaon decay reconstruction method is 
based on the fact that the four dominant K ± decay chan- 
nels (shown in relation |JJ have the same pattern. The 
charged kaon decays into one or two neutral daughter(s) 
which are not detected and one charged daughter which 
is observed in the TPC. 







K 



7T° 



(63.4%) 

TT^+TT (21.1%) 

/*=*= + tt° + ^ (3.27%) 
1 (1.73%) 



(7) 



The decay topology corresponding to the above chan- 
nels is known as a "kink", as the track of the charged 
parent in the TPC appears to have a discontinuity at 
the point of the parent decay. The kink-finding software 
starts by looping over all tracks reconstructed in the TPC 
in the given event, looking for pairs of tracks which are 
compatible with the kink pattern described above. The 
first selection criterion is for the kaon decay vertex (the 
kink) to be found in a fiducial volume in the TPC. The 
TPC has an inner radius of 50 cm and an outer radius 
of 200 cm from the nominal beamline, but the fiducial 
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FIG. 4: Invariant mass distribution of Kg, A, A, H , H + , and D. + Sl + after applying the geometrical cuts outlined in Tables 
□ and|TT| 



volume is defined to have an inner radius of 133 cm and 
an outer radius of 179 cm. The fiducial volume is chosen 
to suppress background due to high track densities (in- 
ner cut) while allowing a reasonable track length for the 
determination of the daughter momentum (outer cut). 
This leads to a maximum number of hits for both the 
parent and daughter track in the fiducial volume. Addi- 
tional cuts are applied to the found track pairs in order 
to select the kink candidates. 

For each kink found, a mass hypothesis is given to both 
the parent and daughter tracks (i.e. K + parent and fi + 
daughter) and the pair invariant mass is calculated based 
on this hypothesis. A cut on the invariant mass (nii nv in 
Table III) can then be applied. As charged pions decay 
with a branching ratio of approximately 100% into the 
same [i + v^ channel as the charged kaons, they will have 
the same track decay topology in the TPC. We therefore 
expect that the kink finding algorithm described above 
will include K + , K~, tt + and 7r~ as kink parent can- 
didates. Therefore, several other cuts must be applied 
to further eliminate the pion background from the kaon 
decays in which we are interested. A summary of the 
applied cuts is given in Table llTTl 



Cut 


K± (kinks) 


Invariant mass 


0.3 < minv < 1.0 GeV/c 2 


Kink angle 


J > asin(p/Ml) + 4. - 1.25 • Ml for p>Ml 
) < asin(p/M2) for p>M2 


Daughter 
mom. 


MOO MeV/c 


DCA/cm 
between 
Parent and 


< 0.123 + 0.082/(p T /GeV/c) 1 - 153 


Daughter 





In Figure El we show the regions excluded by the kink 
angle cut in Table IIIII The cut is placed on the line 



TABLE III: Summary of cuts used in the kink analysis. The 
notation is as follows: 

Ml = (m 2 — m^)/2m M and M2 = (mj — m^)/2m M . See text 
and Figure |K| for further details. 



marking the pion limit. In addition to the cuts listed in 
Table ITTT1 a cut was applied to the specific energy loss to 
remove pion contamination below px=500 MeV/c where 
the kaon and pion dE/dx bands are clearly separated. 

The parent-daughter DCA cut in the last row of Table 
IIII was determined from a two-dimensional (DCA and 
Pt) study of the background. The appropriate cut level 
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FIG. 5: Kink angle cut regions for K + and K identified via 
the kink method. 



was determined in each (DCA,px) cell and the results 
were fit with a function of the form A + (B ■ px )■ The 
resulting parameters A, B, and C are given in Table UTTl 



C. Signal Extraction 

In order to extract the particle yield and (px) , we build 
invariant mass distributions in several px bins for each 
of the particle species except the charged kaons. The 
residual background in each px bin is then subtracted 
through a method referred to here as "bin-counting" . 

In the bin-counting method, three regions are defined 
in the invariant-mass distribution. The first, which is 
defined using the Gaussian signal width found by fit- 
ting the px-integrated invariant-mass distribution with 
a linear function plus a Gaussian, is the region directly 
under the mass peak (±3. 5c, ±4. 5c, and ±2.5<r for the 
Kg, A, and S respectively) which includes both signal 
and background (red or lightly shaded in Figure 0}. For 
the Kg and A invariant mass distributions, the second 
and third regions (blue or dark shading in Figure 0} are 
defined to be the same total width as the signal region 
placed on either side (la away for Kg and A) of the cho- 
sen signal region. For the S, the second and third regions 
are each the size of the signal region and are placed 4cr 
away. In this method the background is implicitly taken 
to be linear under the mass peak. In px bins where the 
background appears to deviate significantly from the lin- 
ear approximation, a second degree polynomial fit is used 
to determine the background under the mass peak. This 
occurs mainly at low px- 

This procedure is carried-out in each transverse mo- 
mentum bin and as a function of event multiplicity. The 
resulting spectrum is then corrected for vertex finding 
efficiency (Section IIII D|) as well as the particle specific 
efficiency and acceptance (Section IIII Afl . The A and A 
spectra are further corrected for higher-mass feed-down 
as detailed in section UlI El 



D. Particle Reconstruction Efficiencies 

The number of reconstructed strange particles is less 
than the actual number produced in the collision due 
to the finite geometrical acceptance of the detector and 
the efficiency of the tracking and decay-finding software. 
Additionally, the quality cuts described in section fill Bl 
reduce not only the combinatorial background but the 
raw signal as well. 

In order to determine the efficiency for each particle 
species as a function of transverse momentum, an embed- 
ding process, similar to that described in section IIII Al 
is employed. In this process, a Monte Carlo generator 
is used to produce the particles of interest with a given 
transverse momentum distribution. The produced parti- 
cles are propagated through the GEANT detector simu- 
lation and the resulting signals embedded into real events 
at the level of the detector response (pixel level). Using 
real events provides a realistic tracking and finding envi- 
ronment for evaluating the performance of the software. 
Only one simulated particle is embedded in any given 
event so as not to overly modify the tracking and finding 
environment. The embedded events are then processed 
with the full reconstruction software chain and the results 
compared with the input to determine the final correction 
factors for the transverse momentum spectra. Whether 
or not the event used for embedding already contained 
one or more strange particles is not a concern as only 
GEANT-tagged tracks are counted for the purpose of 
calculating efficiencies. The resulting total efficiencies 
(acceptance x tracking, finding, and cut efficiencies) are 
plotted in FigureE]for Kg, K ± , A, S, and ft. The correc- 
tion is assumed to be constant over the measured rapidity 
region. 

Finally, a correction needs to be applied to the raw 
particle yields due to low primary vertex efficiencies for 
low multiplicity events described in section If f f Al The 
spectra were binned in multiplicity classes and for each 
class the particle yields were corrected using the proba- 
bilities corresponding to finding a good vertex in an event 
with at least a V0 candidate (open squares in Figure |3J, 
thereby accounting for particles from lost and fake events. 
The overall event normalization is also corrected, using 
the numbers corresponding to the probability of finding 
a vertex (black filled circles in Figure |3J), to account for 
the number of lost events. 



E. Feed-down corrections 

5 and ft baryons produce a A as one of their decay 
products. In some cases, the daughter A can be detected 
as if it were a primary A particle. The result is a modi- 
fication of the measured primary A px spectrum and an 
overestimation of the primary A yield. The amount of 
contamination is unique to the cuts used to find the A. 

In order to correct this, Monte-Carlo S simulations 
were performed and tuned to match the measured shape 
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FIG. 6: Total correction factor (efficiency x acceptance) for Kg, A, K + , K , 3, and £1 after cuts 





Kg 


K± (kinks) 


A(A) 




n + n 


Error Source 


dN/dy 


(Pt) 


dN/dy 


<Pt) 


dN/dy 


<PT> 


dN/dy 


(Pt) 


dN/dy 


(Pt) 


Cuts and Corrections (%) 


5.4 


1.1 


3.7 


2.2 


5.4 


1.3 


13 


1.1 


15 


8.0 


Yield extraction and Fit function (%) 


4.9 


3.7 


1.5 


1.2 


6.3 


4.7 


30 


5.6 


20 a 


3.0 a 


Normalization (%) 


4 


N/A 


4 


N/A 


4 


N/A 


4 


N/A 


4 


N/A 



a The numbers for C+f! are for yield extraction only. The statistics 
do not allow a meaningful fit function study to be done for Q -f CI. 

TABLE IV: A summary of systematic errors from various sources. Errors from yield extraction and fit function for H are from 
comparison between niT-exponential and power-law fits. The normalization error affects only the particle yields. All entries 
are percentages. 
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FIG. 7: Efficiency times acceptance for primary As and As 
coming from H decays. 



and yield of the 2 px spectrum presented in this paper. 
The shape and yield of the A spectrum coming from 5 de- 
cays can then be determined. The total correction factor 



(efficiency x acceptance) was then calculated for both 
primary A baryons and secondary A baryons produced 
by embedded S~ decays (see Figure 0. The correction 
factor is different for A baryons coming from S decays. 
Lastly, the secondary A spectrum is multiplied by the 
correction factor for secondary A baryons, divided by the 
primary A correction factor, and the result is subtracted 
from the measured A spectrum. The application of the 
correction factor is formalized in Equation [SJ 

FD(p T ) = Measured(p T ) - (MC(p T ) ■ R e ff (pt)) (8) 

where FD(pt) is the final feed-down corrected px spec- 
trum, Measured(pT) is the non- feed-down corrected 
Pt spectrum (corrected for efficiency and acceptance), 
A/C(pt) is the secondary px spectrum (determined from 
MC), and i? e //(pT) is the ratio of the secondary efficiency 
and acceptance correction to the primary efficiency and 
acceptance correction. The neutral S° has not been mea- 
sured by our experiment and therefore, for the purposes 
of determining the feed-down correction, the S° yield is 
taken to be equal to the measured S - yield. Similarly, 

the S yield is taken to be equal to the measured S + 
yield. By using Monte Carlo calculations, we determined 
that the finding efficiency for secondary A particles was 
the same whether the A comes from a charged or neutral 
S. Therefore, the final feed-down correction is doubled 
to account for feed-down from S n decays. 
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F. Systematic Errors 

Several sources of systematic errors were identified in 
the analyses. A summary of these errors and their es- 
timated size is to be found in Table IIVI A description 
of various sources of systematic error and their relative 
contribution is given below. 

Cuts and Corrections: The offline cuts that are 
applied to minimize the residual backgrounds also help 
eliminate contamination from pile-up events. The cuts 
may be tightened to further reduce background or loos- 
ened to allow more signal and improved statistics at the 
cost of greater contamination. The final cuts are a com- 
promise between these two extremes which aim to maxi- 
mize the statistics for a given particle species while elim- 
inating as much background as possible. The systematic 
errors from the cut-tuning provide an estimate for our 
sensitivity to changes in the various cuts. 

This number includes the systematic errors from the 
embedding and vertex finding efficiency corrections. The 
A and A entry also accounts for the systematic errors 
from the feed-down correction. 

Methods of Yield Extraction: In order to esti- 
mate the systematic error on the yield extraction in each 
Pt bin, a second method of determining the yield in a 
given bin was used. In the second method, a combination 
of Gaussian plus a linear function is fit to the mass peak 
and background. The yield is then determined by sub- 
tracting the integral of the fitted linear function across 
the width of the signal peak from the sum of the bin 
content in the peak. In both methods, fitting and bin- 
counting, the background is assumed to be linear under 
the mass peak. A second degree polynomial fit is used in 
Pt bins where this assumption is clearly invalid (mostly 
at low pt)- The two methods of extracting the yield may 
give different values due to the finite precision of the fit- 
ting method and fluctuations in the background in the 
bin-counting method. The difference in the two methods 
and any differences resulting from a deviation from the 
linear background assumption are taken into account by 
this systematic error. 

The systematic error on the lowest % 2 /ndf functional 
parameterization of the spectra is also included in this 
row. It is estimated by comparing the yield and 
(pt) from the best % 2 /ndf functional parameterization 
with that from the second best x 2 / ncu " functional form. 
The final numbers for the mid-rapidity yield and (pt) (i n 
Tables fVTI and IVTT|l were determined for each particle us- 
ing the fit with the smallest x 2 /ndf as shown in Table 

El 

Normalization: The overall systematic error from 
the vertex and trigger efficiency affects only the particle 
yields and does not change the shape of the spectra. How- 
ever, the vertex finding efficiency depends on the beam 
luminosity. The number quoted in this row is the level 
of fluctuation in the vertex finding efficiency with beam 
luminosity, 4%. 

Conversion of our measurements to cross-sections must 



also account for an additional 7.3% uncertainty in the 
measured NSD trigger cross-section (26±1.9 mb) and for 
the 86% efficiency of the BBC trigger detectors. 



IV. RESULTS 
A. Spectra 

The fully corrected pt spectra for K + , K~, Kg, A, S, 
and Q are shown in Figure|SJ The measured spectra cover 
only a limited range in transverse momentum and there- 
fore an appropriately parameterized function is needed to 
extrapolate into the unmeasured pt regions for the yield 
determination. In the past, exponential functions such as 
that given in Equation [5] have been used to extrapolate 
spectra from p + p collisions to low transverse momen- 
tum while QCD-inspired power-law functions (see Equa- 
tion [TOJ seem to provide a better description of the high 
p T (>3 GeV/c) region [J Q3, d, EH . The p T coverage of 
the STAR detector for strange particles is large enough 
that a function which accounts for both the power-law 
component of the spectra and the low pt turnover be- 
comes necessary to describe the data. A form that has 
been suggested is the Levy function given by Equation 

\n\M- 



1 



d 2 N 



2npT dydpT 



Ae 



1 d 2 N 
2ttpt dydpr 



d 2 N 



dN 



B{1 + ^) 
Po 



(n-l)(n-2) 



(9) 



(10) 



lirpx dydpT dy 2TmC[nC + mo(n — 2)] 



m 



(11) 



nC 



where A, T, B, p 0l n, ^L, C, and m are fit parame- 
ters. Attempts were made to fit the pt spectra for our 
measured species with all three forms. A summary of the 
resulting x 2 /ndf from each fit is given in Table IVl for each 
of the measured species. The mid-rapidity yields and 
mean transverse momenta quoted below were determined 
from the best fitting form which, for all species, was the 
Levy form (Equation lllfl . The measured mid-rapidity 
yields and feed-down corrected yields are presented in 
Table IVII The measured mean transverse momenta are 
presented in Table IVTT1 

Initially, we compare our measurement of neutral 
strange particles to similar experiments at this energy. 
The closest comparison can be made to the SppS (Super 
Proton- Antiproton Synchrotron) experiments of UA1- 
UA5 using the p + p beam. Only UA5 published strange 
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FIG. 8: Corrected mid-rapidity (\y\ < 0.5) pt spectra for K + , K~, Kg, A, H, and Q. A spectra that have been corrected for 
feed-down are shown as open symbols in the A panel. The dashed lines are fits using Equation II II except for the Q + SI where 
the fit uses Equation |U] The error bars displayed include systematic errors while the fits were done using statistical errors only 
for all species except the charged kaons. 





itit exponential 


Power-Law 


Levy 


Particle 


X 2 /ndf 


ndf 


X 2 /ridf 


ndf 


X 2 /ndf 


ndf 


K° s 


15 


22 


1.5 


21 


0.89 


19 


K+ (kinks) 


3.1 


11 


7.0 


10 


0.40 


9 


K" (kinks) 


9.4 


11 


5.0 


10 


0.30 


9 


A 


4.5 


22 


3.3 


21 


0.81 


18 


A 


4.7 


22 


3.1 


21 


0.99 


18 




0.84 


9 


1.4 


8 


0.76 


8 




1.4 


9 


0.96 


8 


0.83 


8 




0.13 


1 





TABLE V: A summary of \ 2 P er degree of freedom values 
for different fit-functions to the pt spectra. The + fi + 
statistics are considered insufficient for a fit comparison. 



Particle 


dN/dy, \y\ < 0.5 


Stat. Err. 


Sys. Err. 


Kg 


0.134 


0.003 


0.011 


K+ (kinks) 


0.140 


0.006 


0.008 


K~ (kinks) 


0.137 


0.006 


0.007 


A 


0.0436 


0.0008 


0.0040 


A 


0.0398 


0.0008 


0.0037 


A (FD) 


0.0385 


0.0007 


0.0035 


A (FD) 


0.0351 


0.0007 


0.0032 




0.0026 


0.0002 


0.0009 




0.0029 


0.0003 


0.0010 




0.00034 


0.00016 


0.0001 



TABLE VI: A summary of mid-rapidity NSD yields for mea- 
sured strange particles. Numbers in rows marked (FD) have 
been corrected for feed-down as described in section Hll El 



particle measurements at v / s=200 Ge V ll4LIT5ll . with oth- 
ers at Vs=546 GeV [3 and 900 GeV while UA1 at Vs=630 GeV (1] and references cited therein), 
published high statistics strange particle measurements It is worth noting that the UA5 A sample consisted of 
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Particle 


(p T )(GeV/c) 


Stat. Err. 


Sys. Err. 


Kg 


0.605 


0.010 


0.023 


K+ (kinks) 


0.592 


0.071 


0.014 


K" (kinks) 


0.605 


0.072 


0.014 


A 


0.775 


0.014 


0.038 


A 


0.763 


0.014 


0.037 


A (FD) 


0.762 


0.013 


0.037 


A (FD) 


0.750 


0.013 


0.037 




0.924 


0.120 


0.053 




0.881 


0.120 


0.050 


n- +~n + 


1.08 


0.29 


0.09 



TABLE VII: A summary of mid-rapidity (pt) for measured 
strange particles. Feeddown corrected numbers for A and A 
are the same as the non-feed-down corrected values within 
statistical errors. 



only 168 "manually sorted" candidates whereas the 
STAR sample consists of fifty-eight thousand candidates. 

Table lvTlTl compares the values of dN/dy and obtained 
from the STAR pt spectra to the published values from 
the UA5 experiment at SppS |15| measured with a larger 
rapidity interval. In the last column, the STAR data is 
scaled by a factor, obtained via Pythia [ijj simulation, 
to account for the difference in rapidity coverage of the 
two experiments. UA5 measured Kg with |y|<2.5, A with 
|y|<2.0, and 2 with |y|<3.0. STAR measures only in the 
region |y|<0.5. The STAR scaled yields are found to be 
in agreement with the measurement from UA5 and have 
greatly improved on the precision. 

Table Hxl compares the (pt) of the two experiments. 
It was verified, using Pythia, that the dependence of 
(pt) on the different rapidity intervals between STAR 
and UA5 is small, i.e. 2-3%. Therefore, the STAR 
(pt) measurement is compared to UA5 without further 
scaling and is found to have improved on the precision. 



B. Transverse Mass Scaling 

It has been noted previously that the identified particle 
spectra from p+p collisions at ISR energies |2(J,|21| seem 
to sample an approximately universal curve when plot- 
ted versus transverse mass [2^, an effect termed "mx- 
scaling". More recently, data from heavy-ion collisions 
at RHIC have been shown to scale in transverse mass 
over the measured range available |23| . Transverse mass 
spectra from identified hadrons at y / s=540 GeV and 630 
GeV p + p collisions at SppS have also been shown to ex- 
hibit the same behavior up to at least 2.5 GeV j^. The 
degree to which nrr-scaling is applicable and the result- 
ing scaling factors have been used to argue for the pres- 
ence of a gluon-saturated state (color-glass condensate) 
in heavy- ion collisions at RHIC energies [24|. though no 
such interpretation is applied to p + p or p + p collisions. 
Little discussion of the similarity of the results between 



p + p and A+A has been provided. In Figure |9(a)| we 
present the Kg, A, and S mx spectra together with their 
antiparticles and with mi spectra for it, K, and p from 
previously published STAR p + p results at -^=200 GeV 
lHlilEzI- Thc PHENIX 7T° spectrum from p+p colli- 
sions at the same energy is also shown pgj . 

It is clear from Figure 9(a) that while the spectra ap- 



pear to have qualitatively similar shapes, the yields are 
quite different. Nevertheless, the shape similarities en- 
courage us to find a set of scaling factors that would 
bring the spectra onto a single curve. Figure 9(b) shows 
the result of scaling with the set of factors shown in Ta- 
ble [X] These factors were chosen so as to match the it, 
K, and p spectra at an mr of 1 GeV. The higher mass 
spectra are then scaled to match the 7r, K, and p spectra 
in their respective regions of overlap. 

While the low-mx region seems to show reasonable 
agreement between all the measured species the region 
above iit ^2 GeV shows an interesting new effect. The 
meson spectra appear to be harder than the baryon spec- 
tra with as much as an order of magnitude difference de- 
veloping by 4.5 GeV in hit. In order to quantify the 
degree of agreement, a power- law function was fit to all 
the scaled meson and baryon mj spectra separately. The 
ratio of data to each fit was taken for each point in Figure 
9(b)| The data-to-fit ratio is shown for the meson fit in 



Figure [9(c)| and the baryon fit in Figure 9(d) 

This is the first time such a meson-baryon effect has 
been noticed in p + p collisions. This effect is observable 
due to the high px (and therefore high- nix) coverage of 
the strange particle and relativistic rise spectra [271 • The 
harder meson spectrum in the jet-like high-mx region 
may indicate that for a given jet energy, mesons are pro- 
duced with higher transverse momentum than baryons. 
This effect would be a simple reflection of the fact that 
meson production from fragmentation requires only a 
(quark, anti-quark) pair while baryon production requires 
a (di-quark,anti-di-quark) pair. The difference between 
the baryon and meson curves appears to be increasing 
over our measured range, and it will be interesting to 
see, with greater statistics, what level of separation is 
achieved and whether or not the spectra eventually be- 
come parallel. 



Particle Ratios 



Figure 10(a) shows the mean anti-baryon/baryon ra- 
tios (B /B) as a function of strangeness content for p + p 
and Au+Au at yfs = 200 GeV [2<|. The ratios rise 
slightly with increasing strangeness content and are con- 
sistent within errors with those from Au+Au collisions 
at the same center-of-mass energy. Although the B/B 
ratios are not unity for the protons and A baryons, the 
deviation from unity may be explaine d by different par- 
ton distributions for the light quarks 30]. This may be 
sufficient to explain the observed deviation from unity 
without having to invoke baryon number transport over 
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Particle 


STAR dN/dy (|y| < 0.5) 


UA5 Yield 


STAR Yield (scaled to UA5 y) 


Kg 


0.134 ± 0.011 


0.73 ± 0.18, |y| < 3.5 


0.626 ± 0.051 


A + A 


0.0834 ± 0.0056 


N/A 


0.272 ± 0.018 


A + A (FD) 


0.0736 ± 0.0048 


0.27 ± 0.09, jyj < 2.0 


0.240 ± 0.016 




0.0055 ± 0.0014 


0.03±° o: °02, jyj < 3.0 


0.0223 ± 0.0057 



TABLE VIII: A comparison of yields from UA5 (Kg from 15], A from Q] and NSD yields from STAR. The STAR entries in 
the last column have been scaled to the UA5 acceptance using Pythia [l9(l . The STAR errors include systematics. The UA5 
errors shown include their estimated 20% systematic error. 



Particle 


STAR (p T ) (|y| < 0.5) 


UA5 A (p x ) 




0.61 ± 0.02 


0.53t°;^, |y| < 2.5 


A + A 


0.77 ± 0.04 


0.8±°™, |y| <2.0 




0.903 ± 0.13 


0.8i° ;t |y| <3.0 



TABLE IX: A comparison of (p T ) [GeV/c] from UA5 and 
STAR. STAR errors include systematics. The UA5 errors 
shown include their estimated 20% systematic error. 





7T 


K 


P 


A 




Scaling Factor 


1.0 


2.0 


0.6 a 


0.7 


4.0 


Scaled at hit 


1.0 


1.0 


1.0 


1.5 


1.5 



"Data from l25l were scaled by 0.45 

TABLE X: A summary of sca ling fa ctors applied to the trans- 
verse mass spectra in Figure |9(b)| The second row lists the 
transverse mass (in GeV) at which a given particle is scaled 
to match the other spectra. 



five units of rapidity. 

In the case of a quark jet, it is expected that there is a 
leading baryon as opposed to anti-baryon while there is 
no such distinction for a gluon jet. Therefore, making the 
assumption that at high px the observed hadron produc- 
tion mechanisms are dominated by jet fragmentation, it 
is reasonable to expect that the B /B ratio will drop with 
increasing px- This has been predicted previously for cal- 
culatio ns startin g from as low as 2 GeV/c [3l|. Figures 
10(b) and 10(c) show the A/A and 3/3 ratios as a func- 
tion of transverse momentum respectively. Although the 
errors shown in these figures are large, the ratios show no 
sign of decrease in the measured range. The dotted hor- 
izontal line in each figure is the error-weighted average 
over the measured px range. 

One conclusion that could be drawn from the ratios in 
Figure lTUl is that particle production is not predominantly 
the result of quark-jet fragmentation over our measured 
range of px- 



D. Mean Transverse Momentum 

One means of partially characterizing the px spectra 
from p + p collisions is through the determination and 
comparison of the mean transverse momentum. In Figure 
1111 the (px) is shown for all particle species measured in 
both p + p and central Au+Au collisions in STAR. 

In total, twelve particles in both systems are presented, 
covering a mass range of approximately 1.5 GeV/c 2 . The 
solid line is an empirical curve proposed originally 32] to 
describe the ISR and FNAL Q data for it, K , and 
p only, at y/s — 25 GeV. It is interesting that it fits the 
STAR lower mass particles from p + p at y/s = 200 GeV 
remarkably well considering there is nearly an order of 
magnitude difference in collision energy. However, it is 
clear that a different parameterization would be needed 
to describe all of the STAR p + p data. The dependence 
of the inverse slope parameter, T (and therefore of the 
(px)), on particle mass has previously been proposed to 
be due to an increasing contribution to the transverse 
momentum spectra from mini-jet production in p+p and 
p + p collisions j35| . The contribution is expected to be 
even greater for higher mass particles |36|. 

The available statistics allow a detailed study to be 
made. The mid-rapidity px spectra can be binned ac- 
cording to eventwise charged particle multiplicity (un- 
corrected dN c h/dri) and the (px) determined in each bin. 
We present in Figure IT21 the dependence of (px) on un- 
corrected charged particle multiplicity for K + , K , Kg, 
A, and 3. 

The scale difference is readily apparent but perhaps 
more interesting is the increasing trend of (px) with event 
multiplicity. The increase in (px) with multiplicity is 
faster for the A than for the Kg and charged kaons over 
the range from 2 to 6 in (dN c h/dr]}. The statistics avail- 
able in the multiplicity-binned S+S do not allow a proper 
constraint of the Levy fit. The points for S + S shown in 
Figure lT2*l were determined from the error-weighted mean 
of the measured px distribution only. The present level 
of error on the 3 measurement does not allow a strong 
conclusion to be drawn, though the trend with increasing 
mass is suggestive. This mass-ordering of the (px) mul- 
tiplicity dependence has been observed in previous mea- 
surements at three different energies 0. In particular, 
the pions show little increase in (px) when going from 
low to high multiplicity collisions 
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FIG. 9: Comparisons of scaled and un-scaled transverse mass spectra from p + p collisions in STAR and PHENIX at ^=200 
GeV. 7T, K, and p spectra are from |25l l26t |27| while the PHENIX n° spectrum is from |2g|. Error bars include systematics. 



Models inspired by pQCD such as Pythia suggest that 
the number of produced mini-jets (and thereby the event 
multiplicity) is correlated with the hardness (Q 2 ) of the 
collision. The effect of the mini-jets is to increase the 
multiplicity of the events and their fragmentation into 
hadrons will also produce harder px spectra. 



The spectral shape cannot be characterized by a single 
number. It is also possible to compare the multiplicity- 
binncd spectra directly. We show in Figures 13(a) and 
13(b) the ratio (R pp ) of the multiplicity-binned px spec- 



tion ll2fl for Kg and A respectively. 



tra to the multiplicity-integrated (minimum bias) spectra 
scaled by the mean multiplicity for each bin (see Equa- 



-^W(Pt) 



= R 



scale 



> 



(12) 



(13) 



d/V/(ip T (mult, Px) 
<i/V/<ip T (niinbias, p T ) ' 

where 

_ A^ et , e „ t;i (minbias)- < AT e/l (minbias) 
Agents (mult)- < iV c/ \(mult) > 

The changes in incremental shape from one multiplic- 
ity bin to the next then become easier to see. The striking 
change in spectral shape going from the lowest to high- 
est multiplicity bin is further evidence of the increasing 
contribution of hard processes (jets) to the high px part 
of the spectra in high multiplicity events. 
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FIG. 10: Mid-rapidity (\y\ < 0.5) ratios of A to A and S to E~ vs pt- The dashed lines in |10(b)| a nd |10(cJ| are the error- 
weighted means over the measured pt range, 0.882±0.017 for A/A, 0.921±0.062 for S/S. Figure |10(a)| shows the pt averaged 
ra tio for our measured species compared with measurements from Au+Au. The error bars are statistical only. The dashed line 
in |10(a)| is at unity for reference. 
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FIG. 11: (pt) vs particle mass for different particles mea- 
sured by STAR. Error bars include systematic errors. The 
ISR parameterization is given in reference [32| 



MODEL COMPARISONS 



A. Comparison to PYTHIA (LO pQCD) 



Figure 13(c) shows the A/Kg ratio as a function of 
Pt in the various multiplicity bins. We see in all three 
bins that the A shows a sharper increase with px in the 
low px 1.5 GeV/c) part of the spectrum. Furthermore 
there seems to be a relative increase in the A production 
in the intermediate 1.5 < px < 4.0 GeV/c region. 



At the present time, the most ubiquitous model avail- 
able for the description of hadron+hadron collisions is 
the Pythia event generator. Pythia was based on the 
Lund string fragmentation model j33, HH but has been 
refined to include initial and final-state parton showers 
and many more hard processes. Pythia has been shown 
to be successful in the description of collisions of e + e~, 
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FIG. 13: Ratios of multiplicity binned spectra to minimum bias spectra (R PP ) for Kg and A and the ratio of the A spectrum 
to the Kg spectrum in each multiplicity bin. See text for further details. 



T3 

|io- 2 

Jio- 4 
^io- 5 

10 6 
10 7 













■Short 












































• data 






\ PYTHIA 


6.2 




; PYTHIA 


6.3 




: PYTHIA 


6.3, K=3 











0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
P T [GeV/c] 




(a) Pythia for K2. 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
P T [GeV/c] 

(b) Pythia for A. 




(c) Pythia for I 



FIG. 14: Kg (left), A (center), and H (right) pt spectra compared to Pythia (ver 6.22 MSEL1, and ver 6.3) with the default 
K-factor= 1 (solid and dashed curves respectively), and K-factor=3 (dot-dashed curve). 



p + jP and fixed target p + p systems (see for example, ref. 

In this paper we have used Pythia v6.220 and 
v6.317 (using default settings with in-elastic cross-section 
(MSEL=1)) in order to simulate pt spectra for Kg, A and 
S. These have then been compared with the measured 
data. 

As shown in Figure Il4l although there is some agree- 
ment at low pt, there are notable differences above 
PT^l-0 GeV/c, where hard processes begin to dominate. 
Pythia v6.2 underestimates the A yield by almost an 
order of magnitude at pt=3 GeV/c. With the newer 
version 6.3, released in January 2005, these large dis- 
crepancies have been largely reconciled for Kg but re- 
main significant for A and 3. This version includes a 



significantly modified description of the multiple parton 
scattering processes. The red dot-dashed lines in Figure 
El re P r esent a simple tune that was done with Pythia 
6.317 which will be described in more detail below. 

To try and understand the difference between Pythia 
and our results, we made comparisons of (pt) versus un- 
corrected charged multiplicity for Kg and A, as shown in 
Figure IPol As expected from the previous figure, version 
6.2 fails to reproduce the minimum bias magnitude of 
(pt)- Although version 6.3 is capable of reproducing our 
minimum bias values of (pt) it clearly fails to reflect its 
increase with charged multiplicity, suggesting that fur- 
ther tuning is necessary. 

In order to improve the agreement with our data we 
have made some simple changes to the Pythia default 
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FIG. 16: nix-scaling results from Pythia 6.3 with default se ttings. The quark or gluon jet selections are based on the final 
state partons being qq or gg respectively. The second panel ( |16 (b) | also contains mixed final states (qg). The spectra have 
been scaled by the factors listed in the legends. See text for more details. 



parameters. In particular, by increasing the K-factor to 
a value of 3 (set to 1 in the defaults) there is an enhance- 
ment of the particle yield at high px in the model which 
allows it to better describe the data. 

The K-factor, which represents a simple factorization 
of next-to-leading order processes (NLO) in the Pythia 
leading order (LO) calculation, is expected to be between 
1.5-2 for most processes, such as Drell-Yan and heavy 
quark production ^1 & t higher energies. Based on these 
measurements, a K-factor of 3 would signal a large NLO 
contribution, particularly for light quark production at 
RHIC energies. Intriguingly, a large K-factor has been 
estimated for the t/s ~ 200 GeV regime at RHIC based 
on the energy dependence of charged hadron spectra . 
So it seems that for light quark production at lower ener- 
gies, NLO contributions are important and a comparison 



of our data to detailed pQCD based NLO calculations is 
more appropriate. 

With the addition of this K-factor, we can see that the 
Pt spectra for A and 3 in Figure ^] agree even better 
with the model, with the Kg data falling slightly below 
the prediction. More importantly, the Pythia results of 
(px) versus charged multiplicity, including the enhanced 
K-factor, are now in much better agreement with the 
data, as seen in Figure IT31 

Figure El shows the results of separating Pythia 
events based on their final state parton content. Events 
where the final state is qq are labeled as containing 
"quark jets" while events with gg are labled as contain- 
ing gluon jets. Figure 16(a) shows that events with only 



quark-jet final states seem to show a mass splitting in 
the high mx region while events whose final states con- 
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tain jets from gluons (Figure 16(b) I show a shape differ- 
ence between mesons and baryons with the meson spectra 
being harder than the baryon s pectra . The shape differ- 
ence is also apparent in Figure 16 (c)| which contains all 
final states including those with both quark and gluon 
jets. This shape difference could be simply related to the 
fact that a fragmentation process could impart more mo- 
mentum to a produced meson than a produced baryon 
based on mass and energy arguments. This taken to- 
gether with the results shown in Sections If V Ijl and llV CI 
indicates that above 2 GeV in transverse mass, the spec- 
tra contain significant contributions from gluon-)et frag- 
mentation rather than quark-jet fragmentation. 



•k STAR, p+p,\fs.200 GeV 
■ UA1, p+p,>JI=630GeV 

PYTHIA 6.3 default 

PYTHIA 6.3 K=3 




P T GeV/c 

FIG. 17: A/Ks as a function pt compared with Pythia. 



the Kg from Kniehl, Kramer, and Potter (KKP) (43 and 
a calculation by DeFlorian, Stratmann, and Vogelsang 
for the A T4!]. The variations in fi show the theoreti- 
cal uncertainty due to changes of the factorization and 
renormalization scale used. The factorization and renor- 
malization scale allows one to weight the specific hard 
scattering contributions of the parton densities to the 
momentum spectrum. Although for the Kg reasonable 
agreement is achieved between our data and the pQCD 
calculation, the comparison is much less favorable for the 
A. Considering that goo d ag reement was achieved for 
charged pion [23 and 7r° |28l l45j spectra and yields at 
the same energy, our comparison and the comparisons 
in [22) suggest that the region of agreement with NLO 
pQCD calculations may be particle species dependent. 
The baryons are more sensitive to the gluon and non- 
valence quark fragmentation function, which is less con- 
strained at high values of the fractional momentum z 

Recently, the OPAL collaboration released new light 
quark flavor-tagged e + e~ data which allows further con- 
straint of the fragmentation functions |4^ . Albino- 
Knichl-Kramer (AKK) show that these flavor separated 
fragmentation functions can describe our experimental 
data better However, in order to achieve this agree- 
ment, AKK fix the initial gluon to A fragmentation func- 
tion (Dg) to that of the proton (-Dg), and apply an addi- 
tional scaling factor. They then check that this modified 
also works well in describing the p + p SppS data 
at y/s—630 GeV. So, it appears that the STAR data is a 
better constraint for the high z part of the gluon fragmen- 
tation function than the OPAL data. Similar conclusions 
have been drawn elsewhere with respect to the important 
role of RHIC energy p + p collisions 0] . Recent studies 
of forward n° production also suggest that the region of 
agreement with NLO calculations extends as far out as 
3.3 units in (rj) |5Cj . 



In Figure [l3 (c)| we showed the A to Kg ratio separated 
into multiplicity bins. Figure IT7I shows the multiplicity 
integrated ratio compared with Pythia calculations us- 
ing the default settings as well as a K-factor of 3. Here 
we see again the same shape difference between the A 
and the Kg that is seen for baryons and mesons in gen- 



eral in Figure 9(b) and in the p/n ratio 27]. Pythia 
is not able to reproduce the full magnitude of the effect 
in either ratio |27j]. The A to Kg ratio shows a simi- 
lar shape in ^/s = 200 GeV Au+Au collisions though 
the magnitude is larger and multiplicity-dependent |42l| . 
Also, measurements from UA1 at y/s = 630 GeV indicate 
the magnitude may also be dependent on beam energy 

Q 



B. Comparison to NLO pQCD Calculations 

In FigurelTSlwe compare the Kg and A spectra to NLO 
pQCD calculations including fragmentation functions for 



C. Comparison to EPOS 

Finally we compare our data to the version 1.02 of 
the EPOS model |5l| . This model generates the ma- 
jority of intermediate momentum particles by multiple 
parton interactions in the final state rather than frag- 
mentation. The multi-parton cross section is enhanced 
through a space-like parton cascade in the incoming par- 
ton systems. The outgoing, time-like parton emission, 
is allowed to self-interact and to interact with the di- 
quark remnants. The interactions can be either elastic 
or inelastic. The overall result is a strong probability 
for multi-parton interactions before hadronization. The 
cascades are modeled through so-called parton ladders 
which also include multiple scattering contributions of 
the di-quark remnants from a hard parton scattering in 
&p+p collision. Furthermore, by taking into account the 
soft pomeron interactions, the model is able to describe 
the p + p spectra down to low px ■ Finally the inclusion 
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FIG. 18: Kg (left) and A (right) particle spectra (circles) compared to Next-To-Leading-Order (NLO calculations) by Werner 
Vogelsang based on specific Kg |43i| and A (4^| fragmentation functions. Dashed lines illustrate the uncertainty due to the 
choice of factorization scale. 



of parton ladder splitting in asymmetric d+Au collisions 
yields a good description of the difference between p + p 
and d+Au spectra in the same theoretical framework. 
Further details of the model can be found elsewhere [EJ • 

EPOS shows remarkable agreement with BRAHMS, 
PHENIX and STAR data for pion and kaon momentum 
spectra and (px) in p + p and d+Au collisions at both 
central and forward rapidities ( |27l IoTI l52| and references 
therein). Figure ITT?1 shows that this trend also continues 
for the heavier strange particles at mid-rapidity. The 
agreement in p + p collisions in the measured px region 
is largely due to a strong soft component from string 
fragmentation in the parton ladder formalism. Remnant 
and hard fragmentation contributions are almost negli- 
gible at these moderate momenta. The soft contribution 
dominates the kaon spectrum out to 1 GeV/c and the S 
spectrum out to 3 GeV/c. As the momentum differences 
between (di-quark,anti-di-quark) and (quark, anti-quark) 
string splitting are taken into account, and the current 
mass difference between light and strange quarks is folded 
into the spectral shape, a comparison between the spec- 
tra exhibits a flow-like mass dependence. 

The agreement with EPOS is better than even the best 
NLO calculations. A detailed discussion of the differ- 
ences between EPOS and NLO calculations is beyond 
the scope of this paper, but it should be mentioned that 
the two models are, in certain aspects, complementary. 
More measurements of a) heavier particles and b) to 
much higher px are needed in order to distinguish be- 
tween the different production mechanisms. In summary, 
the data show the need for sizeable next-to-leading-order 
contributions or soft multi-parton interactions in order 
to describe strange particle production in p+p collisions. 



D. Statistical Model 

The application of statistical methods to high energy 
hadron-hadron collisions has a long history dating back 



to Hagedorn in the 1960s [53|, 0, EH- Since then sta- 
tistical models have enjoyed much success in fitting data 
from relativistic heavy-ion collisions across a wide range 
of collision energies [5(| |53, IH, Hfl lU |62, |63| ■ The 
resulting parameters are interpreted in a thermodynamic 
sense, allowing a "true" temperature and several chem- 
ical potentials to be ascribed to the system. More re- 
cently, statistical des crip tions have been applied to p + p 
and p+p collisions (64|, and even e + + e~ [6^, but it 
remains unclear as to how such models can successfully 
describe particle production and kinematics in systems of 
small volume and energy density compared to heavy-ion 
collisions. 

It is important to note that a p+p system does not have 
to be thermal on a macroscopic scale to follow statistical 
emission. For example, Bourrely and Soffer have recently 
shown that jet fragmentation can be parametrized with 
statistical distributions for the fragmentation functions 
and parton distribution functions In this picture, 

the apparently statistical nature of particle production 
observed in our data would be a simple reflection of the 
underlying statistical features of fragmentation. It is in- 
teresting to note that Biro and Mueller have shown that 
the folding of partonic power law spectra can produce 
exponential spectral shapes of observed hadrons in the 
intermediate px region with no assumption of tempera- 
ture or thermal equilibrium whatsoever |66|. 

Another possibly related idea is that of phase space 
dominance in which all possible final state configura- 
tions (i.e. those that are consistent with the energy, mo- 
mentum, and quantum numbers of the initial state) are 
populated with equal probability 67]. The finite energy 
available in the collision allows many more final state 
configurations that contain low mass particles than high 
mass particles. The final state configurations containing 
high mass particles are therefore less likely to be observed 
not because they are less probable, but because there are 
fewer of them relative to the low mass configurations. 

We include in this section the results of a canonical 
statistical model fit, using THERMUS jH|, to the STAR 
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FIG. 19: Comparison of Kg , A, H , and 3 ' spectra with calculations from EPOS vl.02. 
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ing THERMUS. Filled circles are ratios from ^=200 GeV 
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tical model fit. All ratios to the left of the vertical line were 
used in the fit. The (Q~ + Q. /2) /S~ ratio was then predicted 
from the fit results. The dashed lines in the lower panel are 
guides for the eye at 1 a. 



feed-down corrected ratios from p+p collisions at y/s=200 
GeV. We used only the canonical formalism as it has been 
determined from a micro-canonical calculation that the 
volume of p+p collision systems does not exceed 100 fm 3 
[69| . Previous results have shown that such a small vol- 
ume invalidates the use of a grand-canonical treatment 
[Tof . The canonical calculation involves only the temper- 
ature (T), baryon number (B), charge (Q), strangeness 
saturation factor (7s), and the radius. For this fit, B and 
Q were both held fixed at 2.0. The resulting parameters 



are presented in Table IXII and a graphical comparison is 
presented in Figure |2*U1 





Canonical Value 


T 


0.1680 ± 0.0081 GeV 


B 


2.000 (fixed) 


Q 


2.000 e (fixed) 


7s 


0.548 ± 0.052 


radius 


3.83 ± 1.15 fm 



TABLE XI: Comparison of a canonical fit to the STAR feed- 
down corrected ratios from p + p collisions at -^5=200 GeV. 
The x 2 /ndf of the fit was 4.14 / 6 = 0.69. See text for further 
details. 



The interpretation of the fit parameters is difficult in 
the context of a p + p collision where the system is not 
expected to thermalize and the volume is small. It is im- 
portant to note that in a pure thermal model, all emitted 
particles would be expected to reflect the same temper- 
ature. Non-thermal effects such as flow would modify 
this result. In p + p collisions, the particle spectra clearly 
show different slopes and those slopes are not in agree- 
ment with the T parameter that results from the statisti- 
cal model fit to the particle ratios. As no flow is thought 
to be present in the p + p system and the results of Sec- 
tion llV Bl supDort that conclusion, this result is a further 
indication of contributions to the particle spectra from 
non-thermal processes like mini-jets. 



VI. SUMMARY AND CONCLUSIONS 



We have presented measurements of K + , K , Kg, A, A, 

S _ , S + , and +f2 + px spectra and mid-rapidity yields 
from i/i=200 GeV p + p collisions in STAR. Corrections 
have been made for detector acceptance and efficiency as 
well as the multiplicity dependence of the primary vertex 
finding and, in the case of the A and A, feed-down from 
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higher mass weak decays. It was found that the measured 
range of transverse momentum necessitates a functional 
form that accounts for the power-law like shape at high 
Pt • We have used a Levy function to fit the spectra and 
extrapolate to low px- 

The (pt) and mid-rapidity yields are in excellent agree- 
ment for all species with previous measurements at the 
same energy but with greatly improved precision. The 
anti-particle to particle ratios are flat with px over the 
measured range for both the A and 5 and therefore show 
no sign of quark-jet dominance at high px- The px inte- 
grated ratios approach unity with increasing strangeness 
content. The anti-baryon to baryon ratios suggest that 
the mid-rapidity region at RHIC is almost baryon-free, 
at least in p + p collisions. The amount of deviation 
from unity expected from differing parton distribution 
functions must first be determined before any claim of 
significant baryon number transport from beam rapidity 
to mid-rapidity can be made. 

We have demonstrated the scaling of transverse mass 
spectra for low px mesons and baryons onto a single 
curve to within 30% out to approximately 1.5 GeV in 
mx- Above 2 GeV the mx spectra show a clear difference 
in shape between mesons and baryons with the mesons 
being harder than the baryons. This is the first obser- 
vation of a difference between baryon and meson spectra 
in p + p collisions and is mainly due to the high px (and 
therefore high mx) coverage of the strange particles pre- 
sented here. Pythia 6.3 seems to account for this effect 
and suggests it is mostly due to the dominance of gluon 
jets. More data are needed to determine the range of the 
effect. 

The mean transverse momentum as a function of par- 
ticle mass from both the p + p and Au+Au systems has 
been compared. Both systems show a strong dependence 
of (px) on particle mass. It is also worth noting that the 
mass-dependence of (px) in the p + p system seems to be 
independent of collision energy as the parameterization 
of the y/s = 25 GeV ISR data seems to work well over 
the same range of measured masses at RHIC. 

The dependence of (px) on event multiplicity was also 
studied for each of the three species (and anti-particles). 
The (px) shows a clear increase with event multiplicity 
for the Kg and A particles. There may be a mass-ordering 
to the increase as the A baryons show a slightly faster 
increase with multiplicity than the Kg, but the present 
level of error on the S data does not allow a definite 
statement to be made. 

The multiplicity-binned Kg and A spectra show a clear 
correlation between high multiplicity events and the high 
px parts of the spectra. The spectral shapes for the 
Kg and A are observed to change with event multiplicity 
and the A to Kg ratio increases over the lower px range 
and reaches higher values in the px range above ~1.5 
GeV/c for larger multiplicites. This suggests that the 
high multiplicity events produce more A hyperons rela- 
tive to Kg than the low multiplicity events. 

Comparisons of our spectra with Pythia v6.221 show 



only poor agreement at best without adjustment of the 
default parameters. In the relatively high px region 
(above 2 GeV/c) there is nearly an order of magnitude 
difference between our data and the model calculation. 
The more recent Pythia 6.3 provides a much better de- 
scription of our Kg data though a K-factor of 3 is required 
to match the A and 5 spectra as well as the observed rate 
of increase of (px) with multiplicity. NLO pQCD calcu- 
lations with varied factorization scales are able to repro- 
duce the high p T shape of our Kg spectrum but not the 
A spectrum. Previous calculations at the same energy 
have been able to match the 7r° spectra almost perfectly, 
which suggests that there may be a mass dependence to 
the level of agreement achievable with pQCD. 

The EPOS model has previously provided excellent 
descriptions of the tt~ , K~ , and proton spectra from 
both p + p and d+Au collisions measured by BRAHMS, 
PHENIX, and STAR at mid-rapidity and forward rapid- 
ity. We have extended the comparison to strange and 
multi-strange mesons and baryons and have found the 
agreement between our data and the EPOS model to be 
at least as good as the best NLO calculations. 

We have demostrated the ability of the statistical 
model to fit our data to a reasonable degree with three 
parameters. Interpretation of the resulting parameters 
in the traditional fashion is not possible as the p + p col- 
liding system is not considered to be thermalizcd. The 
T parameter does not agree with the slopes of the mea- 
sured species and we conclude that this result suggests 
a significant contribution of non-thermal processes (such 
as mini-jets) to the particle spectra. 
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